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Abstract 

We calculate the genuine long-range multi-particle rapidity correlation functions, Cn{yi, ■ ■ ■ ,yn) 
for n = 3,4,5, 6, originating from fluctuations of the fireball longitudinal shape. In these correlation 
functions any contribution from the short-range two-particle correlations, and in general up to 
(n — l)-particle in Cn, is suppressed. The information about the fluctuating fireball shape in 
rapidity is encoded in the cumulants of coefficients of the orthogonal polynomial expansion of 
particle distributions in rapidity. 
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I. INTRODUCTION 


Fluctuations in the longitudinal structure of the hreball produced in heavy-ion collisions 
has drawn noticeable interest in recent years, see, e.g., [1-13]. These fluctuations result in 
new phenomena and modify known correlations in rapidity and azimuthal angle. 

In Ref. [3] it was argued that fluctuations of the hreball longitudinal shape result in the 
specihc long-range two-particle rapidity correlations that depend not only on the rapidity 
difference, yi — 1 / 2 , but also on the rapidity sum, yi + y 2 - In analogy to the long-range 
azimuthal correlations originating from huctuating shape of the hreball in the transverse 
direction [14, 15], it was proposed to expand the single particle rapidity distribution in 
terms of the orthogonal polynomials 


p{y,ao,ai,...) = p{y) 1 + V a^Ti (?/) 

•*—^2 = 0 


( 1 . 1 ) 


where p{y) = {p{y)) is the measured single particle distribution, oq represents rapidity 
independent multiplicity huctuation of the hreball as a whole, oi is an event-by-event asym¬ 
metric component^, etc., see Ref. [3] for more details. Averaging Eq. (1.1) over ao,ai,... 
with the probability distribution P(ao,ai,...) we obtain (oj) = 0. 

The two-particle rapidity distribution at a given oq, ui,... is 


P 2 {yi, 2 / 2 ; ao, ai,...) = p(yi; Oq, Ui, ...)p{y 2 ; ao, Ui,...) -h short range , (1.2) 


where we added the short-range correlations in rapidity, that cannot be written, on an event- 
by-event basis, as a product of huctuating single-particle distributions (1.1). If resonance 
decays or jets are the dominant mechanism behind the short-range contribution, we expect 
this correlation to be approximately a function of yi — y 2 with a range of about one unit in 
rapidity. 

Taking an average over a* and subtracting p{yi)p{y 2 ), we obtain the hnal two-particle 
rapidity correlation function 


^ 2 ( 2 / 1 , 2 / 2 ) = ^ 2 ( 221 , 2 / 2 )-p( 2 /i)p( 2 / 2 ) 

= p(2/i)p(2/2) V., ^{{aiak)+ bik)Ti{yi)Tkiy2) , 


(1.3) 


where bik denotes some short-range (SR) correlations not related to the huctuating shape of 
the hreball 


bik 



Cl^{yi,y2)Ti{yi)Tk{y2) 

p{yi)p{y2) 


(1.4) 


where C'f^(yi,|/ 2 ) is the short-range two-particle correlation function. It is essential to 
remove this unwanted background and this is the subject of the paper. 

We propose to measure the cumulants of the genuine multi-particle rapidity correlation 
functions in analogy to the multi-particle how cumulants [18, 19], which proved to be ehective 
in removing non-how ehects from correlations in azimuthal angle [20]. 


^ In the wounded nucleon model [16, 17] and for symmetric A+A collisions, ai corresponds to the difference 
between left- and right-going wounded nucleons, ai oc wl — wr. 
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Throughout the paper we use y = y, where y is rapidity or pseudorapidity in the 
range [—Y,Y]. The orthogonalization condition for the Legendre polynomials^ [21] reads 

/ dyTi{y)Tk{y) = 6ik, and Tk{y) = ^Jk+\Pk (^) with Pk{x) = - I)''- 

In the next Section we derive formulas for the genuine three-, four-, five- and six-particle 
correlation functions originating from fluctuating longitudinal shape of the hreball. We 
discuss our results in Section 3. 


II. MULTI-PARTICLE CORRELATIONS 

In this section we discuss multi-particle correlations originating from fluctuating hreball 
shape in rapidity. At this point it is useful to comment on the experimental way of estimating 
the integrals over the multi-particle distributions. In the experimental analysis [1], estimates 
of {aittk) are obtained from the integration of the two-particle correlation function. It is 
challenging to apply this method to the higher order cumulants. However, the standard 
procedure of summing over n-tuples, e.g., 

)Tj{yb)Tk{yc)Ti{yd)) = / ^ 2 / 1 '^ 222 ^ 2 / 3 ^ 2/4 - . ^ ^ ^ . 

J-i P(2/i)p(2/2)p(2/3)p(2/4 

' Tjjyg) Tj{yb) Tkjyc) Ti{yd) \ 

ai^,d p^y^'^ p^y<^'^ p^y<^'^ / ’ 

can be applied for samples with sufhcient statistics [22]. P 4 ( 2/1 > 2 / 2 ) 2 / 3 ) 2 / 4 ) is the measured 
four-particle rapidity density. In the last line of the above expression the sum runs over four 
different particles in a given event and the average is over all events.^ 





A. Three-particle correlations 

The three-particle distribution at a given oq, oi,... is 


P3(2/i)2/2)2/3;aO)ai) •••) = p{yi]ao,ai, •••)p( 2 / 2 ; ao) an •••) P(2/3; ao) an •••) + short range. (2.2) 


Expanding on the orthogonal basis and taking an average over a* we obtain 


P3(2/1)2/2)2/3) 

p(2/i)p(2/2)p(2/3) 


1 + ((fliOfc) + hk) (2/1) Tk (2/2) + Ti (2/1) Tk (2/3) + Ti (2/2) Tk (2/3)] + 

{{aiakam) + bikm)Ti{yi)Tk{y2)Tm{y3) , (2.3) 


where bik and bikm represent the unwanted short-range background (not related to fluctu¬ 
ations of the fireball longitudinal structure). We are interested in extracting information 


^ In Ref. [3] we expanded the distribution (1.1) in the Chebyshev polynomials but other choices are certainly 
possible. 

^ By definition {Ti{ya)) = f dyTi{y) ~ 
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about the genuine three-particle correlations'^, C 3 (|/i, 2 / 2 , 2 / 3 ), dehned as 

P3(2/1, ?/ 2 , 2 / 3 ) = p(2/i)p(2/2)p(2/3) + p( 2 /i)C 2 ( 2 / 2 , 2 / 3 ) + p( 2 / 2 )C' 2 ( 2 /i, 2 / 3 ) + 

p(2/3)C'2(2/i, 2 / 2 ) + <^ 3 ( 221 , 2 / 2 , 2 / 3 ), (2.4) 


where ^3(2/1,2/2,2/3) is the three-particle rapidity density. If genuine three-particle short-range 
correlations can be neglected, = 0, performing simple calculations we obtain 


<^ 3 ( 2 / 1 , 222 , 2 / 3 ) = p(2/i)p(2/2)p(2/3) 




i,k,m=0 


Ti (^l) T]^ (2/2) (l/s) 


( 2 . 5 ) 


Equation (2.5) allows to directly extract {aittkam) without the contribution from the 
two-particle short-range correlations bik 


/ dyidy2dy^ 


<^ 3 ( 2 / 1 , 2 / 2 , y^)T^{yl)Tk{y 2 )Tm{y■i) 
p(2/i)p(2/2)p(2/3) 


( 2 . 6 ) 


where (...)[ 3 ] denotes^ that {aittkam) is sensitive to but it does not depend on the lower 
order correlation function 6 * 2 . 

Using Eq. (2.4) we can relate ( 0 * 0 ^ 0 ™) [3] through integrals of the two- and three-particle 
densities pn and hnally through sums over n-tuples, see Eq. (2.1) 


{ttittkttrr) j-gj 


{Ti{ya)Tk{yb)Tm{yc)) - {Ti{ya)) {Tk{ya)Tm{yb)) - (TkiVa)) {Ti{ya)Tm{yb)) - 
{Tm{ya)) {Ti{ya)Tk{yb)) + 2 {Ti{ya)) {Tk{ya)) {Tm{ya)) 

(Ei(2/a)Tfc(2/6)r,n(2/c))[3] . ( 2 . 7 ) 


Perhaps U 3 is not particularly useful because (af) is expected to be rather small if not 
zero (by dehnition {a) = 0). Specihc effects of correlations resulting from bias of event 
multiplicity on rapidity distribution, e.g., stronger longitudinal expansion in events with 
higher hreball density, can be tested using mixed correlations 

(aoafc)[3] = {To{ya)Tk{yb)Tkiyc)) - {Toiya)) {Tk{ya)Tk{yb )), (2.8) 

where {T^iya)) = \/2 for the chosen normalization and P > 0 . 


B. Four-particle correlation 


Here we discuss the more interesting case of four-particle correlation function. Performing 
analogous calculations we obtain 

= 1 + W. ^ „ {a,a,) [T, („0 n (!/ 2 ) + T, (»,) T, (y,) + D (Si) U {v,} + 
P(2/i)p(2/2)p(2/3)p(2/4) ^^,k=o 


Ti (2/2) Tk (2/3) + Ti (2/2) Tk (2/4) + Ti (2/3) Tk (2/4)] + 
y] {aittkttm) [Ti (2/1) Tk (2/2) Tm (2/3) + Ti (2/1) Tk (2/2) Tm (2/4) + 

^ ^2,fe,m=U 

Ti (2/1) Tk (2/3) Tm (2/4) + Ti (2/2) Tk (2/3) Tm (2/4)] + 
y (aiafcOman) r* (2/1) (2/2) Thi (2/3) T'n (2/4) . ( 2 . 9 ) 

^ ^i,k,m,n=0 


^ In general, the genuine n-particle correlation function C„(?/i, ■■■,yn) is non-zero only if there is a physical 
mechanism directly correlating n or more particles. 

^ For C 3 we have (aiCfcOmlp] = {a-iO-kam), which is not the case for higher order correlation functions. 
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For simplicity we omit the short-range background correlations bik, hkm and bikmn- We 
are interested in the genuine four-particle correlation function, (^ 4 ( 2 / 1 , ?/ 2 , ?/ 3 , 2 / 4 ), dehned as 


2/2,2/3,2/4) = p(2/i)p(2/2)d(2/3)d(2/4) + p(2/i)p( 2/2)C'2(2/3, 2/4) + p(2/i)p(2/3)C'2(2/2,2/4) + 
p(2/i)p(2/4)(^2(2/2,2/3) + p(2/2)p(2/3)(^2(2/i, 2/4) + p(2/2)p(2/4)C'2(2/i, 2/3) + 
p(2/3)p(2/4)(^2(2/i, 2/2) + p(2/i)C'3(2/2, 2 /3, 2/4) + p(2/2) <^3(2/1,2 /3, 2/4) + 
p(2/3)C'3(2/i, 2 /2, 2/4) + p(2 /4)(^3(2/1,2/2,2/3) + (^2(2/1,2 /2)C'2(2/3, 2/4) + 
C'2(2/l,2/3)C'2(2/2,2/4) + (^2(2/1,2/4)<^2(2/2,2/3) + (^4(2/1,2 /2,2/3, 2/4)- (2.10) 


Performing straightforward calculations we obtain 


<^ 4 ( 2 / 1 ,2/2, 2/3, 2 / 4 ) 
p(2/i)p(2/2)p(2/3)p(2/4) 


E (ninfcnmnn) [41 ( 2 / 1 ) Ffc ( 2 / 2 ) ( 2 / 3 ) ( 2 / 4 ), 

2,fc,m,n=0 ^ J 


where 


( 2 . 11 ) 


[4j — iflk^rn) • 


( 2 . 12 ) 


with (...)[ 4 ] denoting that the object depends on (74 but not on the lower order correlations. 

At this stage we are mostly interested in extracting the leading terms, i = k = m = n 
(and in particular the asymmetric term oi) 

^ {at)-S{a^y 


<^4(2/1, 2 / 2 , 2 / 3 , 2 / 4 )Ti (2/1) Ti (2/2)7; (1/3) Ti (1/4) 

= / dyidy2dy^dyi - r \ r \ r \ r \ - 

J P(2/i)p(2/2)p(2/3)p(2/4) 

= {Tiiya)Ti{yb)Ti{yc)Ti{yd)) - 3{Ti{ya)Ti{yb))‘^ , 


(2.13) 


where in the last line of the above equation (/ > 0) we show the way to calculate the 
cumulant, see Eq. (2.1). 

Another interesting term in the expansion (2.11) is the mixed term (ao^l) k > 0. The 
expression for such a cumulant reads 

(aga^)[4] = {alal) - (al) {al) - 2 {aoakf 

= {To{ya)To{yb)Tk{yc)Tkiyd)) - {To{ya)To{yb)){Tk{ya)Tk{yb)) - 
2{To{ya)Tk{yb)f - 2{To{ya)){To{ya)Tk{yb)Tk{yc)) + 
2{To{ya)f{Tk{ya)Tk{yb)). (2.14) 


This expression removes two and three-particle correlations and correctly takes into account 
that (To(2/a)) 7 ^ 0. In particular (aQa|)[ 4 ] could be a measure of the genuine correlations 
between event multiplicity and the width of the particle distribution in rapidity. 


C. Five-particle correlation 


For the genuine hve-particle correlation function, < 75 ( 2 / 1 , ..., 2 / 5 ), dehned as® 
p5 = PPPPP + pCi + PpC '3 -|- pppC2 + PC 2 C 2 + <72173 -|- < 75 , 


( 2 , 16 ) 


10 


10 


15 


10 


For clarity we skip the argument yi and show only the numbers of possible combinations. 
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we obtain 


<^5(yi, •■•, 2 / 5 ) 

p{yi)-p{yh) 

where 


rgi Ti (?/l) Tfi (?/2) (2/3) (2/4) (Z/s) ) 

2,/c,m,n,r=0 ^ 

(^2) [5] — iflm^n^r) “t“ •■•]• 

10 variations 


( 2 . 16 ) 

( 2 . 17 ) 


The leading term is 

(af)[5] ^ {4)-l^{aj){a\) (2.18) 

= {Ti{ya)Ti{y^)Ti{y^)Ti{yd)Ti{y^)) - lQ{Ti{ya)Ti{yb)){Ti{ya)Ti{y^)Ti{yc)), 


where in the second line of the above equation we assume i > 0. 


D. Six-particle correlation 


Finally, for the six-particle correlation function, C^iih, ...,2/6), dehned as 

Pe = PPPPPP + P^b + PpCi + PPpCs + PPPpC 2 + PC2C2, + + C2C4 + 

6 15 20 15 60 45 15 

C^+C^^+Ce, (2.19) 

10 15 


we obtain 

^* 6 (^ 1 , y&) 


E,., 


CLiQj}^(lm^n^r^s) [g] Ti (t/l) Tj^ (2/2) (^/s) Tji (2/4) T^ (Z/s) Tg (l/o); 

( 2 . 20 ) 


p(^l)...p(^6) ^i,k,m,n,r,s=0 

where 


15 variations 


(^n^r^s) ~l~ •••] ~l~ (^r^s) "h •••]• (2.21) 

^ - 


10 variations 


The leading term is 


a^)[6] = (a^) - 15(a^)(af) - lO(a^)^ + 30(a)^)^ 


3\2 


15 variations 


.2\3 


( 2 . 22 ) 


III. COMMENTS AND CONCLUSIONS 

We propose to measure higher order cumulants of event-by-event fluctuations of rapidity 
distribution. The multi-particle rapidity distributions p„(i/i,..., yn) can be expanded in the 
basis of orthogonal polynomials Ti^{yi) ... Ti^{yn) with coefficients (a^j ... ai^}. The coeffi¬ 
cients have contributions from the rapidity shape fiuctuations of the distribution function. 
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that we are after, and from the short-range correlations. Calculating higher cumulants of 
such averages reduces the contribution from short-range correlations. For example, in the 
fourth order cumulant {Ti{ya)Tj{yb)Tk{yc)Ti{yd))[ 4 ] the contribution from two- and three- 
particle short-range correlations is removed and, if the fourth-order short-range correlations 
are neglected, it can be directly compared to the fourth cumulant of the expansion coeffi¬ 
cients (ajajafca;)[ 4 ]. 

The extracted value of the higher cumulants ... ai„)[n] can be compared to predictions 
of models of event-by-event fluctuations of rapidity distributions (1.1). The cumulant of the 
coefficients a* can be calculated in models once their event-by-event distribution is known. 
The proposed method to study the cumulants of the expansion coefficients does not rely on 
the precise model assumptions for fluctuations of rapidity distributions. It remains a subject 
of further studies to calculate all significant four or six order cumulants, in different models 
of energy deposition in hadronic collisions. The direct calculation of higher order cumulants 
from experiment or Monte Carlo events requires a very larger statistics. The limited statistics 
available has prevented us from applying the method to realistic hydrodynamic events for 
Pb-Pb collisions at the LHC. 

Finally, we would like to emphasize that our method is applicable not only to symmetric 
A-l-A collisions but also to asymmetric interactions including p-|-A. It would be also inter¬ 
esting to perform measurement in p-|-p collisions, where the internal quark (diquark [23]) 
structure of a proton should results in, e.g., nonzero asymmetric term (aj‘)[„], n = 2,4,6. 

Acknowledgments: 

Research supported by the Ministry of Science and Higher Education (MNiSW), by found¬ 
ing from the Foundation for Polish Science, and by the National Science Centre (Narodowe 
Centrum Nauki), Grant Nos. DEC-2014/15/B/ST2/00175, DEC-2012/05/B/ST2/02528, 
and in part by DEC-2013/09/B/ST2/00497. 


[1] The ATLAS collaboration [ATLAS Collaboration], ATLAS-CONF-2015-020; 
http://cds.cern.ch / record /2029370 

[2] A. Bialas, A. Bzdak and K. Zalewski, Phys. Lett. B 710, 332 (2012) [arXiv: 1107.1215 [hep- 
ph]]. 

[3] A. Bzdak and D. Teaney, Phys. Rev. C 87, no. 2, 024906 (2013) [arXiv:1210.1965 [nucl-th]]. 

[4] P. Bozek, W. Broniowski and J. Moreira, Phys. Rev. C 83, 034911 (2011) [arXiv:1011.3354 
[nucl-th]]. 

[5] P. Bozek, W. Broniowski and A. Olszewski, Phys. Rev. C 91 , 054912 (2015) [arXiv: 1503.07425 
[nucl-th]]. 

[6] A. Bialas and K. Zalewski, Nucl. Phys. A 860, 56 (2011) [arXiv:1101.1907 [hep-ph]]. 

[7] A. Olszewski and W. Broniowski, Phys. Rev. C 92 , no. 2, 024913 (2015) [arXiv: 1502.05215 
[nucl-th]]. 


7 



[8] Y. Cheng, Y. -L. Yan, D. -M. Zhou, X. Cai, B. -H. Sa and L. P. Csernai, Phys. Rev. C 84, 
034911 (2011) [arXiv:1106.3371 [hep-ph]]. 

[9] L. P. Csernai, G. Eyyubova and V. K. Magas, Phys. Rev. C 86, 024912 (2012) [arXiv: 1204.5885 
[hep-ph]]. 

[10] L. Pang, Q. Wang and X. -N. Wang, Phys. Rev. C 86, 024911 (2012) [arXiv; 1205.5019 [nucl- 
th]]. 

[11] V. Vovchenko, D. Anchishkin and L. P. Csernai, Phys. Rev. C 88, no. 1, 014901 (2013) 
[arXiv:1306.5208 [nucl-th]]. 

[12] J. Jia and P. Huo, Phys. Rev. C 90, no. 3, 034915 (2014) [arXiv; 1403.6077 [nucl-th]]. 

[13] L. G. Pang, G. Y. Qin, V. Roy, X. N. Wang and G. L. Ma, Phys. Rev. G 91, no. 4, 044904 
(2015) [arXiv; 1410.8690 [nucl-th]]. 

[14] S. A. Voloshin, A. M. Poskanzer and R. Snellings, arXiv;0809.2949 [nucl-ex]. 

[15] B. Alver and G. Roland, Phys. Rev. C 81 , 054905 (2010) [Erratum-ibid. C 82 , 039903 (2010)] 
[arXiv;1003.0194 [nucl-th]]. 

[16] A. Bialas, M. Bleszynski and W. Czyz, Nucl. Phys. B 111 , 461 (1976). 

[17] A. Bialas and W. Czyz, Acta Phys. Polon. B 36, 905 (2005) [hep-ph/0410265]. 

[18] N. Borghini, P. M. Dinh and J. Y. Ollitrault, Phys. Rev. C 63, 054906 (2001) [nucl- 
th/0007063]. 

[19] N. Borghini, P. M. Dinh and J. Y. Ollitrault, Phys. Rev. C 64, 054901 (2001) [nucl- 
th/0105040]. 

[20] V. Khachatryan et al. [CMS Collaboration], Phys. Rev. Lett. 115 , no. 1, 012301 (2015) 
[arXiv;1502.05382 [nucl-ex]]. 

[21] J. Jia, S. Radhakrishnan and M. Zhou, arXiv; 1506.03496 [nucl-th]. 

[22] A. Bilandzic, R. Snellings and S. Voloshin, Phys. Rev. C 83 , 044913 (2011) [arXiv; 1010.0233 
[nucl-ex]]. 

[23] A. Bialas and A. Bzdak, Phys. Rev. C 77, 034908 (2008) [arXiv;0707.3720 [hep-ph]]; Phys. 
Lett. B 649, 263 (2007) [nucl-th/0611021]. 


